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SUMMARY 


An analytical study has been made of the effects of several airplane and 
landing-gear design variables on the airplane response to runway roughness (for a 
relatively smooth runway) with application to supersonic transport configurations. 
Assumptions made for the analysis were (l) natural frequency of all tires were 
the same, (2) tires had no damping and possessed a linear tire -deflection load 
curve, (3) all struts had linear damping and spring characteristics, and the 
natural frequency and damping of all struts were the same, (^) rigid airplane 
structure, and ( 5 ) no aerodynamic forces. 

The results of the analytical study indicated that, within the limitations 
of the assumptions made, the parameter variations considered had little effect on 
motions of the center of gravity of the airplane except for longitudinal position 
of the landing gear with respect to the center of gravity; in the latter case, 
moving the landing gear aft with respect to the center of gravity caused a reduc- 
tion in center-of-gravity displacement and acceleration at all speeds. The 
effects of parameter changes on motions of the pilot's compartment were strongly 
influenced by speed; that is, changes which increased the response at one speed 
could cause decreased response at another speed. From the standpoint of root- 
mean-square accelerations at the pilot's compartment, placing the main gear near 
the center of gravity appeared to be advantageous; however, it resulted in 'some- 
what higher displacements than with the main gear placed further rearward from 
the center of gravity. Generally, the supersonic transport configurations gave 
higher root -mean-square accelerations at the pilot's compartment than did a sub- 
sonic jet transport configuration. 


INTRODUCTION 


As presently envisaged, supersonic transports are expected to have an unusu- 
ally long fuselage, a major portion of it extending forward of the center of 
gravity. This design will allow the possibility of landing-gear arrangements 
and locations relative to the pilot and passenger compartments substantially 
different from those of current aircraft. During take-offs and landings, these 


configurations may cause motions or accelerations undesirable for cockpit instru- 
ments, pilots, and/or passengers. 

The present study was undertaken to examine the effects of such differences 
on the responses of the pilot's compartment and the airplane center of gravity 
due to runway roughness of a relatively smooth runway. The results are also com- 
pared with estimates for a current subsonic jet transport airplane. 


SYMBOLS 


a,b,c 


C 

C S 

D 

D* = D/a 
F 


g 



k s 


k t 

k y 

i 


*n 

L 


m 

M 

s 
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constants 

constant 

strut damping coefficient 
differential operator, d/dt 

input force from runway 

acceleration due to gravity, 32.2 ft/sec^ 

strut spring constant 
tire spring constant 
pitch radius of gyration, ft 

wheel base (nose gear to main gear), ft 
distance from center of gravity to main gear, ft 

distance from center of gravity to nose gear, ft 

wavelength, ft 

static mass on one gear 

static mass on all wheels 

distance along runway, ft 

airplane speed along runway, ft /sec 

longitudinal distance from center of gravity, positive forward, ft 
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generalized output term 


generalized input term 


vertical displacement, ft 
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1 + — 

K, 


8 ^ J ky2 J ky 2 




_ J*cg 


K 2 v V 


pitch angle, radians 

damping ratio of strut with static load mass 


root-mean- square displacement, ft 


root -mean- square acceleration, g 


root -mean- square pitch angle, deg 


phase angle, radians 

power- spectral-density function, sq ft/radian/ft 

power-spectral-density function, g^/radi an/ft 
frequency, radians/sec 
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natural frequency of strut with static load mass, 20 radians/sec 


Subscripts : 


natural frequency of tire with static load mass, 10 radians/sec 

natural frequency of oscillation in pitch 

frequency ratio, (jo/a>t 

spatial frequency, 2 jc/L, radians/ft 

center of gravity 


input 

refers to point between tire and strut 


main gear 


nose gear 


associated with oscillation in pitch 
output 

pilot compartment 


runway 


strut 


x at x-position 

Dots over symbols represent differentiation with respect to time. Bars ovei 
symbols denote vector quantities . 


ASSUMPTIONS, ANALYSIS, AND METHOD OF EVALUATION 


The following assumptions were made: (l) the natural frequency of the tires 

with the mass equivalent to the static load, was the same for the main and nose 
wheels, (2) the tires had no damping and possessed a linear tire-deflection- load 
curve, (3) all struts had linear damping and spring characteristics and the natu- 
ral frequency and the damping of all the struts were the same, (4) the airplane 
fuselage was a rigid structure, and (5) there were no aerodynamic forces. 

The equations of motion of the airplane resulting from forces applied to the 
tires as the airplane moved along the runway were derived under the assumptions 


k 


stated. The resulting frequency response was combined with the power spectrum of 

6.7 x 10~ 6 

power spectrum of acceleration at various points along the airplane longitudinal 
axis resulting from center-of- gravity translations and pitching responses. The 
root -mean-square accelerations, displacements, and pitch angles were determined 
by integrating these spectra over the range of wavelengths from 4 feet to 
570 feet. The computations were carried out by using a digital computer for 
50 points in this range for each of the 7 discrete speeds of 40, 80, 120, 160, 

200, 240, and 280 feet per second. A complete development of the equations used 
in the analysis is presented in the appendix. Figure 1 defines the dimensional 
symbols representing the characteristics of the airplane . 

The parameters which were varied were landing-gear wheel base Z, longitudi- 
nal location of the landing gear with respect to the center of gravity, and radius 
of gyration ky . Responses were calculated for various points along the longi- 
tudinal axis. 

The long slender fuselage envisaged for the supersonic transport may be more 
flexible than those of current jet transports. The acceleration responses for 
such a fuselage may be somewhat different from those presented herein for an 
assumed rigid fuselage, depending on such factors as the modes excited, their 
natural frequencies, arrangement of the landing gear, airplane speed, and so forth. 


(see ref. l) to give the 


runway roughness for a good runway $(fi) r 


Presentation of Results 

The results of the analytical study are presented in figures 2 to 4. Table I 
shows the airplane conf igurations , values of the parameters used, and the figure 
numbers in which these various configurations appear. Figure 2 shows the varia- 
tion with velocity of the root -me an- square values of normal accelerations , the 
root -me an -square normal displacements, and the root -me an- square pitch angles due 
to the variation of the landing-gear location with respect to the center of gravity 
(fig. 2(a)), the variation of the pitching radius of gyration (fig. 2(b)), and the 
variation of the length of the landing-gear wheel base (fig. 2(c)). Figure 3 shows 
the variations of root-mean-square accelerations and displacements with location 
along the longitudinal axis for various speeds for a possible supersonic transport 
configuration having a short wheel base (45 ft) relative to overall length 
(fig. 3(a) ) and for a 90-foot wheel base (fig- 3(h) ) . Figure 4 compares the vari- 
ation with velocity of the same quantities shown in figure 2 for the supersonic 
transport and a current subsonic turbojet transport calculated by the same method. 


DISCUSSION 


Effect of Landing-Gear Location 

The effect of varying the longitudinal location of the landing gear with 
respect to the center of gravity so that the main gear varied from 0.1Z to 0.32 
aft of the center of gravity (fig. 2(a)) is to reduce the root-mean-square dis- 
placements and accelerations of the center of gravity at all speeds. The 
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root-mean-square acceleration at the pilot ' s compartment nearly doubled at the 
higher speeds as the landing gear was moved aft; however, the maximum accelera- 
tion was only about 0.3g and this was for a relatively smooth runway. From the 
standpoint of root-mean- square accelerations at the pilot's compartment, placing 
the main gear near the center of gravity appeared to be advantageous; however, 
it resulted in somewhat higher displacements ( 0.3-foot maximum) than with the 
main gear placed further rearward from the center of gravity (about 0.2-foot 
maximum displacement). The root-mean- square pitch angles were less than 0.2° for 
all these gear locations and showed no important variations with speed. 


Effect of Pitching Radius of Gyration 

As would be expected, changing the radius of gyration in pitch had no effect 
on the values of root-mean- square accelerations and displacements at the center of 
gravity. Halving the pitching radius of gyration from 28.3 feet (ky/Z = O.63) to 
14.15 feet (ky/Z = 0.315) resulted in a threefold increase in root -mean- square 
accelerations at the pilot's compartment at the hi^ier speeds, with the maximum 
value about 0.45g. (See fig. 2(b).) 

At a speed of 240 feet per second it may be noted that, for the pilot's com- 
partment, although the values of the displacements and the pitch angles are about 
equal for the various configurations, there is a substantial variation in the 
values of root -mean- square accelerations. The reason, at least in part, is prob- 
ably the difference in the natural frequency in pitch of the configurations. The 
values of root -mean- square accelerations should be proportional to the square of 
the natural frequencies since the displacements are essentially the same. 


Effect of Length of Landing-Gear Wheel Base 

Increasing the length of the landing-gear wheel base (by moving the nose 
wheel forward) while also maintaining a constant pitch radius of gyration showed 
that substantial reductions could be realized in both root-mean-square pitching 
angles and root-mean- square displacements for the pilot's compartment over almost 
the entire range of velocities. (See fig. 2(c).) For example, doubling the 
wheel base from 45 feet to $K) feet results in a reduction of root-mean-square dis- 
placements by a factor of as much as 6 at a velocity of about 120 feet per second 
and a reduction of root -me an- square pitch angle by a factor of about 5 for the 
same speed. Increasing the landing-gear wheel base had mixed effects on accelera- 
tions, depending on speed. Although the intermediate wheel-base length (67.5 feet) 
yielded values of root-mean- square pitch angles and displacements which fell gen- 
erally between those for the 45- and 90-foot wheel bases, this intermediate wheel- 
base length produced the greatest root -mean- square accelerations at speeds above 
120 feet per second. 


Variation With Distance Along the Fuselage 

For the configuration with a wheel base of 45 feet (fig. 3(a) )> the values 
of displacement and acceleration response at speeds of 40, 80, 120, and 160 feet 
per second appear to be a mi nimum at or near the center of gravity and exhibit a 
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fairly smooth and uniform increase in the values with increasing distance either 
toward the nose or tail at all speeds. For the configuration with a wheel base 
of 90 feet, however (fig. 3(h) ) , the minimum values of displacements and acceler- 
ations occur at points somewhat displaced from the center of gravity (x/i = 0) 
for all speeds (in this figure, 40, 80, 160, and 2^0 ft/sec). For a speed of 
kO feet per second the minimum acceleration was near the center of gravity; at 
80 feet per second the minimum response was toward the tail (x/i = - 0 A) ; and for 
speeds of 160 and 2b0 feet per second the minimum response locations were toward 
the nose (x/i ** 0 * 5 ). 


Comparison of Hypothetical Supersonic Transport Configuration 
With Current Subsonic Jet Transport 

The curves of figure 4 indicate that, for the pilot's compartment for the 
^ 5 -foot wheel base, the root-mean-square displacements of the supersonic transport 
were greater than those for the subsonic transport by factors from 1 ^ to 2 ^ over 

the speed range. However, for the 90-foot wheel-base supersonic transport, the 
displacements were about the same or somewhat lower than those for the subsonic 
transport. The root -mean- square acceleration and pitch angles for the supersonic 
transport (^ 5 -foot wheel base) were up to 50 percent larger than those for the 
subsonic transport. For the 90-foot wheel-base configuration, the pitch angles 
were lower than those for the subsonic transport over the entire speed range; 
however, the root-mean-square accelerations were higher than those for the sub- 
sonic transport at the lowest and highest speeds and about the same in the mid- 
speed range (110 to 180 ft/sec). 


CONCLUDING REMARKS 


A simplified analysis has been made to examine certain design variables of 
landing- gear location and aircraft pitch radius of gyration in relation to pos- 
sible effects on the response of supersonic transport configurations to runway 
roughness . The results indicate that the parameter variations considered had 
little effect on motions of the center of gravity of the airplane except for 
longitudinal position of the landing gear with respect to the center of gravity; 
moving the landing gear aft resulted in a reduction in center-of -gravity dis- 
placement and acceleration at all speeds. The effects of the parameter changes 
on the motions of the pilot's compartment were strongly influenced by speed; that 
is, changes which increased the response at one speed could cause decreased 
response at another speed. From the standpoint of root-mean- square accelerations 
at the pilot ' s compartment, placing the main gear near the center of gravity 
appeared to be advantageous; however, it resulted in somewhat higher displace- 
ments than with the main gear placed further rearward from the center of gravity. 

Increasing the pitching radius of gyration tended to reduce the root-mean- 
square accelerations and increase displacements at the pilot's compartment. 
Increasing the landing-gear wheel base tended to decrease pilot's compartment 
displacement response but had mixed effects on accelerations, depending on speed. 
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The highest accelerations were obtained with an intermediate wheel base of 
67.5 feet in a range from feet to 90 feet at speeds above 120 feet per second. 

Generally, the supersonic transport configurations gave higher root-mean- 
square accelerations at the pilot's compartment than did a subsonic jet transport 
configuration . 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 2, 1962. 
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APPENDIX 


THEORETICAL ANALYSIS 


Determination of Acceleration Response of the Airplane 
Due to Runway Roughness 

The following assumptions are considered in this analysis. For each wheel, 
the static load is considered to he rated static load, and for this condition the 
natural frequency on all tires is assumed to he the same . In addition, the 

strut response is assumed to he linear, and all struts have the same frequency 
responses and damping ratios (with static load on given wheel as a separate mass 
and neglecting mass of wheels and struts) . 


Determination of the Transfer Function for the Center-of -Gravity 
Translation and Pitching of the Airplane 
With reference to figure 1 the following equations may he written: 


F = k t ( Zi - Zl ) = k s ( Zl - z 0 ) + C s ( Zl - z 0 ) 


(1) 


The equations of motion in vertical translation are: 


MS 


! 'cg 0 ' F n + F m = *t n (zi - *l) n + *^(*1 - sj 


m 


( 2 ) 


Since 


“n = 


1^4 and 
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-> 


ie C g 

L °o m n 1 v x 


z-.) + — — - Zq) 

i/n mm l \ 1 ^m 


(3) 


1m. Isi 

If — 5. = — -3. = then 


m r 


% 


^cg Q = <at 


! {V( Z i ' + t( Z 1 ' z l) 


m 


(4) 
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Also 


m cg 0 ' k s n (*l - z o)„ + Cs n (il - io) n + *^(*1 - ^o) m + c s m (il - *o) m (5) 


or 


Mz 




S-n Zr 


c §r 


+ _ z j + _^21 ^Wz-, - Zq) 

m n i \ x °J m nu 1 V 1 °/ 


m 


(6) 


If -^1 = — = (D g ^ and -^S. - = 2Agcn s , then 

m n 3% %i 3% 


Z =8o ’ "s 


2 ^*(*1 - *o) n + 2V“s ^(*1 - io) n 


+ m s 2 - z 0 ) m + 2^% ^(il - i 0 ) 


m 


(7) 


The following relation also may he written: 


c Si,!,o = 7^1,1,° + TV,l,o 


( 8 ) 


Substituting equation (8) into equations (4) and (j) yields 


z cg 0 = ^ 2 ( z i - z ]) 


eg 


(9) 


and 


= m_2 


( Z 1 - z o) c + 2A s m s (z 1 - z 0 ) ( 


Z nrr = (1)„ 

c &o s \ -*■ eg ° ° \ -L w / eg 


( 10 ) 
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The equation of motion in pitch is 


Mky 9 q — x n F n - (11) 

Substituting the expressions for F n and F m from equation (2) into equa- 
tion (ll) , together with the equations for 1 %, anfl cDj. 2 , yields 


Mky^e'o = M^2o>t 2 ( Z;L - Zi) n - M^Sto fc 2 (z i - zj 


( 12 ) 


From figure 1, = ~ — and Therefore, 


Also 


® 0 - ‘ 9 l) 

*y 


(13) 


Mky^e'o = M^£o) s 2 (z 1 - z D ) n - (zq - z Q ) m 


+ 2M^HVo s (z 1 - z Q ) n " (zi - z 0 ) m 


(l*) 


or 


So - ^ 0 2 K - e o ) + 2^FV>s(9l - «o) 

k y V 


(15) 


The solution of equations (9), ( 10 ) , (13), and (15) yields the transfer functions 
for the center -of- gravity translation and pitching responses of the airplane to 
runway inputs. Equations ( 9 ) and ( 13 ) are of the form 


X D = a 2 (x 1 - Xj 

and equations (10) and ( 15 ) are of the form 


( 16 ) 


*o = * 2 (*1 - *o) + c(*q - X G ) 


(17) 
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where a, h, and c are constants, and X^ is the displacement of the runway 
undulations or the input, and X Q is the resulting displacement of the airplane 
or the output. is the displacement relating to motions of the wheel axles 

and lower parts of the shock struts. Writing equations (l6) and (17) in operator 
form and rearranging results in the following: 

D 2 X q + a 2 X-L = a 2 X i (l8) 


(d 2 + cD + b 2 )x o - (cD + b 2 )^ = 0 (3 

Equations (l8) and (19) can he evaluated simultaneously hy eliminating X^_ to 
obtain 


D 3 + fa 2 ,±b 2 V + a 2 D + a^b£ 


X n = a 


l%> + ^^|X 


Taking D = aD* and rewriting equation (20) yields the nondimensional form: 


5.1+ (b/a) n* 2 j. n* 


D + 


+ V* + 


X„ = D* + 


: / a P 


where 


+ aD* 2 + D* + ^X Q = (d*'+ p)x i 


1 + (b/a) 2 


p = ( W *) 2 

c/a 


From equation (22) the response ratio is obtained: 


- oui* 2 ) + cd* 2 ^ - m* 2 ) 2 + (a - 0)<n* 5 
^0 - arn * 2 ^ 2 + a .* 2 ^1 - ca * 2 ^ 2 
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where 


or 


and the phase angle between output and input: 




Z Cg ± 



■where 


and 


“eg = “/“t 


a cg - 


1 + kh ) 2 

2^s(“s/“t) 


Peg 


_ “W®t 


2A 0 


(a,p,(o*) 0 


(a,p,co*) cg 



The next step is to determine the input information z cg ^ and 0^ in terms of 

the runway characteristics. If it is assumed that the runway roughness spectrum 
is made up of sinusoidal waves of spatial frequencies C and the runway dis- 
placement at the nose wheel is taken as reference, the input motion at the nose 
wheel for unit amplitude is 


Z *i 



( 27 ) 


where 


s is distance along the runway and the input at the main wheels is 





ift(s-Z) 


(28) 


The displacement at the center of gravity resulting from displacements at the 
nose and main wheels is 
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z 


cgi 


nu 

= IV 



or 



Is + Ing-iflJ ] e ifts 

l l J 


from which 


J cg ± _ ^c gi _ 


II - £Vn(i _ c 


Jn i 


Z2 


and the phase angle is 


n sin ill 

0 ogl = sin' 1 1 


2Z m Z, 


II - -^Kl - CO 


The pitch -angle displacement input is 


e i^ - z n i “ z m i 


so that, from equations (27) and (28), 


hi » (i - e -isn) e i» s 

z n ± ' ' 


from which 


= J2(l - cos S1Z) 

* y * 


z n^ z r 


0 fl = sin-1 sin 


2(1 - cos ftZ) 


and the phase angle is 



or 


^9± “ “ til) 


(55) 


Equations ( 27 ) to (35) can be transformed to the time domain by the relation 

a = cs/v. 

The phase angles of the aircraft motion relative to the runway displacement 
can then be determined from equations ( 2 6 ), ( 31 ) , and (35) as 


^°cg,9 (^° " ^ eg, 0 + ^ i cg,0 


(56) 


The motion of the airplane at any distance x from the center of gravity 
along the x-axis where x is positive forward is then given by 


% _ Zg gQ , X l9 P 


l TT 


(57) 


or 


Co . ^#<*0 + x V *>< 
- z r * z_ 


ifl *0e 


( 58 ) 


The response ratio at x is then given by 


If ■ \Wf + 1 tj ^ (w) 


The acceleration response at any point along the x-axis is then 


z T z r 


(ho) 


The power spectral density of airplane vertical acceleration in response to 
runway roughness is 

t(a) s - (lr) * (fl)r 


m 
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where the runway roughness spectrum was taken according to reference 1 to he 


®(a) r 



(te> 


with the value of C = 6.7 X 10 ^ representing a good runway, used for the 
computations . 

The root-mean-square acceleration d* z - at various points along the x-axis 
of the airplane were obtained from the relation 

2 * 
r Lf 

°2 2 - 4t “ 

Lo 


2rtV 

Lf 


.1 r 

V J 2rtV 


dCD 


For the speeds considered in the analysis and considering the frequency- 
response characteristics of the airplane there appeared to he little power in 
the airplane-acceleration response spectrum for runway wavelengths greater than 
570 feet; L 0 was therefore given this value. For the upper limit of the inte- 
gration, the wavelength Lp was taken as k- feet because runway roughness 
measurements do not ordinarily go below this value. 
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TABLE I.- VALUES OF PARAMETERS USED IN ANALYSIS 


If ft 

D 

D 

ky/l 


radians/ sec 

Figure 

Supersonic configuration 

^5 

0.1 

0.9 

0.63 

28.3 

4.28 

2(a) , 2(b), 2(c), 4 

^5 

.2 

.8 

.63 

28.3 

5.72 

2(a) 

^5 

.3 

.7 

.63 

28.3 

6.62 

2 (a) 


.1 

•9 

.50 

22.5 

5.40 

2(b), 3(a) 

45 

.1 

.9 

.315 

14.15 

8.76 

2 (b) 

^5 

.1 

•9 

•50 

22.5 

5 . 4 o 

3 (a) 

^5 

.1 

•9 

•50 

22.5 

5 . 4 o 

3 (a) 

45 

.1 

•9 

.50 

22.5 

5.40 

3 (a) 

45 

.1 

•9 

.50 

22.5 

5.40 

3 (a) 

67.5 

.1 

•9 

.419 

28.3 

6.50 

2(c) 

90 

.1 

•9 

.515 

28.3 

8.76 

3 (b) 

90 

.1 

•9 

.515 

28.3 

8.76 

3 (b) 

90 

.1 

.9 

•315 

28.3 

8.76 

3 (b) 

90 

.1 

•9 

•515 

28.3 

8.76 

2(c), 3(b) 

90 

.1 

•9 

.515 

28.3 

8.76 

3 (b) 

Subsonic jet transport 

52.53 

.1 

.9 

•50 

26.16 

5.40 

4 
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;ure 5.- Variation of root-mean- square normal acceleration and no rmal displacement 
with location along the longitudinal axis of a supersonic airplane configuration 
for various speeds. l m /Z = 0.1; l n /l = 0 . 9 . 









